Mon. Not. R. Astron. Soc. 000, 000-000 (2005) Printed 5 February 2008 (MN M?eX style file vl.4) 



An empirical model for the polarization of pulsar radio 
emission 

Don Melrose, Andrew Miller, Aris Karastergiou^ & Qinghuan Luo 

School of Physics, University of Sydney, NSW 2006, Australia 
— Received in original form 

ABSTRACT 

We present an empirical model for single pulses of radio emission from pulsars based 
on gaussian probability distributions for relevant variables. The radiation at a specific 
pulse phase is represented as the superposition of radiation in two (approximately) 
orthogonally polarized modes (OPMs) from one or more subsources in the emission 
region of the pulsar. For each subsource, the polarization states are drawn randomly 
from statistical distributions, with the mean and the variance on the Poincare sphere 
as free parameters. The intensity of one 0PM is chosen from a log-normal distribution, 
and the intensity of the other 0PM is assumed to be partially correlated, with the 
degree of correlation also chosen from a gaussian distribution. The model is used to 
construct simulated data described in the same format as real data: distributions of 
the polarization of pulses on the Poincare sphere and histograms of the intensity and 
other parameters. We concentrate on the interpretation of data for specific phases of 
PSR B0329-I-54 for which the OPMs are not orthogonal, with one well defined and 
the other spread out around an annulus on the Poincare sphere at some phases. The 
results support the assumption that the radiation emerges in two OPMs with closely 
correlated intensities, and that in a statistical fraction of pulses one 0PM is invisible. 
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1 INTRODUCTION 

r Recent observations of the polarization of single radio pulses 
■ from pulsars (Karastergiou et al. 2001, 2002; Karastergiou 
, et al. 2002; Karastergiou, Johnston & Kramer 2003) have 
> confirmed much earlier observations (Ekers & Moffat 1968; 
Clark & Smith 1969; Lyne, Smith & Graham 1971) that the 
polarization can vary substantially from pulse to pulse and 
that the radiation appears to be a mixture of two orthogo- 
nally polarized modes (OPMs) (Stinebring et al. 1984; McK- 
innon & Stinebring 2000; McKinnon 2002; McKinnon 2003a; 
Johnston 2004) . The emphasis in the recent observations has 
been on the polarization at a specific pulse phase (Karaster- 
giou, Johnston & Kramer 2002; Karastergiou, Johnston & 
Kramer 2003; Edwards & Stappers 2004; Edwards 2004; 
Johnston 2004) , whereas the emphasis in the early literature 
was on the variation with time within a single pulse. Ob- 
servations of the Stokes parameters, I,Q,U,V, at selected 
pulse phases provide statistical information on four parame- 
ters at each phase, and these are conventionally represented 
by the intensity, /, the degree of polarization, r, and a point 
on the Poincare sphere, describing the state of polarization, 
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represented by the latitude, 2x, and longitude, 2i/). For each 
chosen phase, the polarization in each pulse is represented 
by a point on the Poincare sphere, and the OPMs corre- 
spond to a distribution of points concentrated around the 
ends of a diagonal through the sphere. Several different pro- 
jections of the Poincare sphere have been used to present 
pulsar data: a Mercator projection (Lyne, Smith & Graham 
1971), a Hammer- Aitoff projection (Karastergiou, Johnston 
& Kramer 2003), and a Lambert equal area projection of 
the hemispheres (Edwards & Stappers 2004; Edwards 2004). 
The representation of the data on the Poincare sphere con- 
tains information on only two of the four Stokes parameters, 
and to be complete it needs to be complemented with infor- 
mation on the statistics of the intensity and the degree of 
polarization, and on the correlations between these two pa- 
rameters and the polarization. Such information is available, 
and is usually represented by histograms. 

The pulse to pulse variation in polarization is probably 
an important clue to understanding the emission process and 
its polarization characteristics. The interpretation in terms 
of OPMs is very strongly indicative of an emission process 
that generates radiation in two orthogonal modes with corre- 
lated intensities. The polarization appears to vary randomly 
about a mean for each OPM, with the ratio of the intensi- 
ties also varying randomly. In order to identify the physical 
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processes that lead to these variations it is important to 
describe them statistically. McKinnon & Stinebring (2000) 
proposed a model in which the polarization of each OPM is 
always 100%, and the observed polarization is < 100% due 
to it being the sum of contributions from the two OPMs. 
Johnston (2004) fitted data to this model, assuming gaus- 
sian variations with a mean and a variance for the relevant 
variable. Concentrating on the circular polarization, McKin- 
non (2002) proposed a statistical model, that includes noise 
convolved with the signal, and found that large fluctuations 
in circular polarization are possible due to the variations in 
relative intensity of the two modes. McKinnon (2003a) de- 
veloped this statistical model further by including the other 
Stokes parameters in it, leading to a probability distribution 
for the polarization on the Poincare sphere, and McKinnon 
(2004) applied this formalism to observed fluctuations in 
the position angle (PA) of the linear polarization. Edwards 
& Stappers (2004) also used a model that involves a prob- 
ability distribution on the Poincare sphere in interpreting 
their data on PSR B0329-I-54. Our model incorporates all 
these features with the exception of noise, and adds features 
relating to nonorthogonality of OPMs and the correlation 
between the intensities of the OPMs. 

In the empirical model described here, we assume gaus- 
sian probability distributions to describe each relevant vari- 
able in the model. The basic model involves the emission 
being the sum of the two OPMs each of which has its polar- 
ization described by a two-dimensional probability distribu- 
tion on the Poincare sphere. This involves four parameters 
for each OPM: two to describe the mean (e.g., the mean 
latitude and longitude on the Poincare sphere) and two to 
describe the variance, which we assume to be different in two 
different directions. There is observational evidence that the 
intensity variations from pulse to pulse are consistent with 
log-normal statistics (Cairns, Johnston & Das 2001), with 
notable exceptions such as giant pulses. The statistics are 
log-normal at different phases, although the statistical pa- 
rameters can vary with pulse phase (Cairns, Johnston & 
Das 2004). In the model we specify that the intensity of one 
OPM, identified as mode 1, has log-normal statistics. The 
mean of In 7i is set to zero, without loss of generality. The 
probability distribution then involves one additional param- 
eter, the variance in In Ji, whose value, Ai = 0.8, we choose 
to match the observed variance in the natural log of the 
intensity. The interpretation of the data requires that the 
intensities of the two OPMs be neither strictly correlated or 
completely uncorrelated: they must be partially correlated. 
We model the partial correlation by assuming that the in- 
tensity in mode 2 is fc times that in mode 1 (in a given pulse) 
with k described by a gaussian distribution with a mean fco 
and a variance A2. A notable result is that we find that some 
features require fco close to but not equal to unity, and A2 
small, implying a tight correlation. 

Noise is not included in our gaussian probability distri- 
butions; it was included in probability distributions by McK- 
innon (2002; 2003a; 2004). We could include noise by con- 
volving each relevant probability distribution with a noise 
distribution. However, this would increase the number of free 
parameters in the model, and complicate the identification 
of the important parameters needed to explain the features 
on which we concentrate here. Our neglect of noise implies 




Figure 1. The polarization of single pulses of PSR B0329+54 
(at 328 MHz) are shown for two of eight specific phases in Fig. 2 
of Edwards & Stappers (2004). Each panel contains the Lambert 
projection of each hemisphere of the Poincare sphere, with the 
colour code showing the density of points (one point per pulse). 
All eight lead to significantly different distributions, with only 
the first and last (shown here in the lower panel) having obvious 
simple interpretations in terms of OPMs. The example shown in 
the upper panel has one reasonably well defined mode (left hemi- 
sphere) with the other mode seemingly spread out along an an- 
nulus on the sphere. The lines in the lower panel are for constant 
latitude, 2x, and longitude, 2-!/;. 



that the model applies only to data for which instrumental 
noise and other noise is known to be unimportant. 

To illustrate the use of the model we focus on the in- 
terpretation of the data presented by Edwards & Stappers 
(2004) for eight different phase bins of PSR B0329-I-54. For 
some of these bins the data are well described, at least to a 
first approximation, by an OPM model with slightly ellip- 
tically polarized modes and some random variations in the 
polarization of each OPM, with the means being orthogonal. 
For other phase bins the OPMs are not orthogonal, and in 
several examples one mode is spread out around an annu- 
lus on the Poincare sphere. We focus on the data that seem 
most difficult to explain within the model: nonorthogonality 
of the OPMs, and asymmetry between the OPMs. 

The probability distributions are introduced in Sec- 
tion 2, and the results of specific simulations are presented in 
Section 3. Physical interpretations are discussed in Section 
4, and the conclusions are given in Section 5. 



2 PROBABILITY DISTRIBUTIONS 

The model is based on gaussian probability distributions for 
each of the relevant variables: these distributions are defined 
in this section. 
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2.1 Representation on the Poincare sphere 

A state of polarization is represented by a point on the 
Poincare sphere. The longitude, 2^p, on the sphere deter- 
mines the position angle of linear polarization, PA = tp. 
Similarly, the latitude, 2x, determines the axial ratio of the 
polarization ellipse, T = tanx: right hand circular, T — 1, 
corresponds to the north pole, 2x = 7r/2, and left circular, 
T = —1, to the south pole, 2x = —n/2. Points in between 
the pole and the equator represent elliptical polarizations. 
In terms of the Stokes parameters, it is convenient to write 
q = Q/I, u = U/I, V — V/I, so that the degree of polariza- 
tion is r = (g^ + + v^Y^'^ . It is only the polarized part 
that can be represented by a point on the Poincare sphere, 
and this corresponds to (Melrose & McPhedron 1991) 

- = cos(2x) cos{2tp), — = cos(2x) sm{2tp), 



= sin(2x). 



(1) 



Information on I and r is not included in this representation, 
and must be presented in some complementary way. 

A specific example of the representation of the polar- 
ization of individual pulses on the Poincare sphere is shown 
in Figure 1, and is discussed further in Section 3. The num- 
ber of pulses with polarization close to a given point on the 
sphere may be represented as a probability density for the 
polarization represented by that particular point. 

2.2 Polarization probabilities 

In the model the emission in a single pulse is assumed to be 
a sum of completely polarized components, consisting of two 
OPMs for each of N subsources. The polarization state for 
a particular 0PM and particular subsource is drawn from 
a two-dimensional gaussian distribution on the Poincare 
sphere. This formalizes a suggestion made by Karastergiou, 
Johnston & Kramer (2003) that, from pulse to pulse, the po- 
larization vectors vary randomly about a mean. The mean 
polarization is represented by a point, X = = V'i say, 

with i = 1,2 labeling the two modes. The mean polarizations 
may also bo described by qi,Ui,Vi, which are related to Xi^^i 
by (1) with r = 1. We assume that the probability distribu- 
tions around these means may be treated independently for 
the two modes. The probability distribution around mode i 
may be plotted on the Poincare sphere in terms of contours 
of constant probability. The data are suggestive of distribu- 
tions that are anisotropic, favoring one direction over an- 
other, and to allow for this we assume that the contours of 
constant probability arc ellipses. If these surfaces were circles 
we would need only one parameter to describe the variance 
in the polarization for mode i, and with the generalization to 
ellipses of constant probability we require three parameters 
to describe the variance for each mode (the extra two being 
say the axial ratio and the direction of the major axis of the 
ellipse). To minimize the number of parameters, while still 
allowing for auisotropy, wc adopt the following prescription 
that fixes the orientation of these probability ellipses. 

Let the axes on the Poincare sphere be rotated such 
that the new north pole coincides with the mean polarization 
point for mode i, and let quantities relative to the rotated 
axes be denoted by primes. Thus this choice corresponds to 
v'i = 1, q'i = u'i = Q. li there were no anisotropy, a gaussian 



probability distribution with variance a\ would be propor- 
tional to exp[— (g'^ -I- u''^)/2a^\. Wo include an anisotropy 
and fix its orientation by choosing probability distributions 
in q' and u' with different variances, af and 6?, respectively. 
With this prescription, the q' are drawn randomly from the 
distribution 



-g'^/2a2 



-1< g' < 1 



(27raf)V2' 
and the u' are drawn randomly from 

Pu' = 7:^-rT;TT777, -Ku < 1. 



(27r62)i/2' 
Then one has 
V = {1 — q — u ) ' . 



(2) 



(3) 



(4) 



The two-dimensional distribution on the sphere is then rep- 
resented by expressing the product of the probabilities (2) 
and (3). The ratio of the semi-axes of the ellipses of constant 
probability is equal to ai/hi. On rotating back to the original 
axes, this determines the orientation of probability ellipses 
in terms of the original (unprimed) variables centered on 
the mean polarization state (2xi, 2il)i). For linearly polar- 
ized modes {qi = ±1, Ui = Vi = 0), these ellipses are oriented 
along latitude and longitude on the Poincare sphere. In the 
cases discussed below the mean polarization is nearly linear, 
and these probability ellipses are then oriented nearly, but 
not exactly, along latitude and longitude. 



2.3 Nonorthogonal OPMs 

Strictly orthogonal OPMs correspond to antipodean points 
on the Poincare sphere: their PAs differ by 90° and their 
degrees of circular polarization are equal in magnitude and 
opposite in sign. This corresponds to 2x2 ~ —2xi, 2ip2 = 
2%p-i -\- TT. Nonorthogonal OPMs correspond to nonzero dif- 
ferences 2Ax = 2(x2 + Xi), SAi/) = 2{^2 - ipi) - n. These 
two differences are free parameters in the model. 



2.4 Intensities of the two modes 

One mode, mode 1 say, is assumed to have a log-normal 
intensity distribution: 



Ph = 



-(ln/i-lni-o)V2A^ 

(27rA2)i/2 



(5) 



where In Jo and Ai are the mean and the variance of In Ji , 
respectively. The intensity of mode 2 is chosen to be pro- 
portional to the intensity of mode 1, with the constant of 
proportionality, k, selected randomly from a gaussian distri- 
bution with mean ko and variance A2: 



kh, 



-(fe-fco)V2A^ 



Pk 



(2^Ai; 



1/2 



(6) 



This ensures that the intensities of the two modes are par- 
tially correlated, as suggested by observations (Stinebring et 
al. 1984; Karastergiou, Johnston & Kramer 2002). 
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2.5 Number of subsources 

One interpretation of the microstructure in pulses is that 
at any given instant the observer received radiation from 
a number of separate subsources (Luo 2004). Also, using 
a carousel model to explain drifting subsources, Edwards, 
Stappers & van Leeuwen (2003) argued that some data sug- 
gest that two subsources might be visible simultaneously. 
Our empirical model allows us to regard each pulse as com- 
posed of emission from A'^ separate subsources, with different 
choices of parameters for each subsource. To avoid introduc- 
ing too many free parameters, in the simulations reported 
here, when there is more than one subsource, all the sub- 
sources have the same probability distributions. Each sub- 
source then corresponds to an independent set of random 
numbers chosen according to the specified probability dis- 
tributions. Even in this simple case, the results for N sub- 
sources are not related as simply to that for a single sub- 
source as might be anticipated. 



2.6 Additional parameters 

A specific difficulty arises in explaining phases where one 
of the OPMs is concentrated around a well-defined mean 
polarization and the other is spread out around an annulus 
on the Poincare sphere. An additional assumption seems to 
be required to account for these two features being present 
simultaneously. We allow this by introducing an additional 
parameter which is a probability that mode 2 is present 
and mode 1 is invisible. A physical effect that allows this is 
that the ray for mode 2 from a given subsource reaches the 
observer, but the ray for mode 1 from that subsource does 
not because it is refracted into a different direction (Allen & 
Melrose 1982). For each pulse, the model randomly decides 
whether to add the orthogonal modes from all the subsources 
or only a fraction of them. This adds a new parameter to 
the model: Pv, the fraction of pulses for which mode 1 is 
invisible. For multiple subsources, N > 1, we assume that 
Pv applies only to a fraction, pat, of the subsources. 



3 RESULTS OF SIMULATIONS 

We present six examples of our simulations. The parame- 
ters chosen for these six cases are summarized in Table 1. 
Others parameters not listed, P — 15000 (the number of 
pulses in each simulation), Iq ~ 1.0, Ai = 0.8, 2'i/>o = 90°, 
2xo ~ —5° and ai = 0.05 are the same for all the sim- 
ulations reported here. In presenting the results of specific 
simulations, we start with the case of a single subsource with 
orthogonal OPMs, then relax the assumption of orthogonal- 
ity and the assumption of a single subsource. To focus our 
simulations we attempt to identify parameters such that our 
results simulate the specific observations shown in Figure 1. 



3.1 Method of calculation 

In attempting to find a model that simulates the data shown 
in Figure 1, we start by making the simplest assumptions of a 
single subsource (A'^ = 1) with orthogonal OPMs (2Ai/) — 0, 
2 Ax ~ 0). With this first example, the polarizations of the 
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Figure 2. An example of the output of the model is shown. The 
top panel shows the polarization on the Poincare sphere in the 
same format as Figure 1. The polarization points in the right 
hemisphere are dominated by mode 1 and those in the left hemi- 
sphere by mode 2. The first three histograms show the intensity 
summed over the set of simulated pulses: the first panel shows the 
total intensity (in the color version the total intensity is in black 
and its components in two OPMs are in red and green), and the 
second and third panels show the linearly and circularly polar- 
ized intensities, respectively. The fourth panel shows the PA of 
the linear polarization, and the final two panels show the degrees 
of linear and circular polarization. In this example the OPMs are 
orthogonal and their variances are the same and independent of 
direction on the Poincare sphere. There is only one subsource. 
The ratio of the intensity of mode 2 to mode 1 has a mean of 
unity and a variance of 0.5. The other specific parameters are 
given in the text. 



two modes are orthogonal, with equal variances that are in- 
dependent of orientation on the Poincare spheres, and with 
the same mean intensity. The simulated data are shown in 
Figure 2. The first panel is a Lambert projection of the 
two hemispheres of the Poincare sphere, in the same format 
as used in Figure 1. The hemispheres are roughly centered 
on the two mean polarizations, mode 1 in the right hemi- 
sphere and mode 2 in the left hemisphere. The information 
on the intensity is shown in the first of the six histograms. 
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Table 1. The parameters used in the simulations shown in Figures 2 to 7 are tabulated. Angles are in 
degrees. The parameters P = 15000, lo = 1.0, Ai = 0.8, 2^)0 = 90°, 2x0 = -5° and ai = 0.05 are the same 
in all figures, and the parameter ppf is not applicable (NA) for a single subsource, TV = 1. 
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Figure 3. Similar to Figure 2 with only the Poincare sphere 
shown; the centroids of the OPMs remain orthogonal but the 
variances are different in different directions on the sphere and 
different for the two OPMs (details in the text). 



with the intensities in the two modes shown separately in 
red (mode 1) and green (mode 2). The second and third 
and fourth histograms show the linearly polarized intensity 
(L = {Q^ + U^)^^'^), the circularly polarized intensity and 
the distribution in position angle. The usefulness of these 
histograms is for comparison with data, which can be rep- 
resented as histograms of I,L,V and PA with minimal pro- 
cessing. The final two panels show the degree of linear and 
circular polarizations. We construct such histograms for all 
our examples, but here we show only the Poincare sphere 
for most of them. 

The simulation shown in Figure 2 reproduces the exam- 
ple shown in the lower panel in Figure 1 reasonably well. It 
is clear that this and similar looking examples may be simu- 
lated in terms of two orthogonal modes with random spreads 
about the two polarizations. However, it is also clear that 
additional assumptions need to be made to simulate the ex- 
ample shown in the upper panel in Figure 1. We attempt 
to do so by relaxing one-by-one some of the assumptions 
made in the simulation shown in Figure 2. Also note that 
despite the variances in the polarizations being the same, 
and the mean intensities being the same, for both modes, 
the distributions on the Poincare sphere are different; this 
difference is particularly notable in the histogram for PA. 
These differences are due entirely to the different statistics 
for the intensities of the two modes: in this case the highest 
intensities are dominated by mode 2, with {I2) = (Ii) and 
(/|> = 1.25(/?>. 



3.2 Orthogonal OPMs, different variances 

The example illustrated in Figure 3 differs from that in Fig- 
ure 2 in that the variances in the polarizations of the two 
modes are different, and the mean of the ratio of the in- 
tensities of the two modes is increased to slightly greater 
than unity (fco — 1.1). For mode 1, the choice ai = 0.05, 
bi = 0.15, implies that the spread is larger on the Poincare 
sphere in the vertical than the horizontal direction; the sur- 
faces of constant probability are ellipses with axial ratio 
bi/ai — 3, favoring circular over linear polarization. For 
mode 2, the choice 02 = 0.02, 62 = 0.01, corresponds to a 
smaller spread than for mode 1, and such that the proba- 
bility ellipses have axial ratio fe2/a2 = 1/2, favoring linear 
over circular polarization. The effect, shown in Figure 3, is 
to spread out the points around the mean for mode 1 (on 
the right) much more strongly than for mode 2 (on the left). 
It is clear that while such spreading might be one ingredi- 
ent in attempting to simulate the upper panel in Figure 1, 
it cannot explain the concentration of the points around a 
broad annulus, rather than a central modal point. The in- 
crease in fco from unity in Figure 2 to fco = 1.1 in Figure 3 is 
relatively unimportant; this ratio becomes more important 
in Figures 4 to 7 where we attempt to simulate the annulus 
in Figure 1. 

The example illustrated in Figure 4 differs from that 
in Figure 3 in that the spread in the ratio of intensities is 
much smaller: the parameter A2 — 0.5 in Figure 3 is re- 
placed by A2 = 0.05 in Figure 4. This causes the distri- 
bution of points to become even more strongly spread out, 
and for the concentration of points around the means to 
disappear. In interpreting this, first consider a case (not 
shown) where there is no spread in the ratio of the in- 
tensities of mode 2 to that in mode 1, fco = 1, A2 — > 0. 
Then, on summing over the Stokes parameters for the two 
modes the mean polarizations cancel. The cancellation is 
not exact because the parameters for each mode correspond 
to two different choices of random numbers. The net po- 
larization is determined by the difference between the two 
modes, and the degree of polarization is necessarily small, 
[{Qi + Q2f + iUi + U2f + {Vi + V2fV^^ < I1+I2 = 2/1, due 
to {Q2, f/2, V2) ~ —Qi, Ui,Vi. This leads to a broad distri- 
bution of polarization on the Poincare sphere. The difference 
between Figure 3, A2 = 0.5 and Figure 4, with A2 = 0.05 is 
that the latter is effectively indistinguishable from the case 
A2 ^ in which the mean polarizations cancel exactly. 

This example adds a further ingredient that is plausi- 
bly needed in the interpretation of the upper panel in Fig- 
ure 1: spreading out of the points due to near equality of 
the intensities in the two modes. However, the associated 
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Figure 4. As for Figure 3 but with the variance in the ratio 
of the intensity of mode 2 to mode 1 reduced from A2 = 0.5 
to A2 = 0.05 so that there is a strong tendency for the mean 
polarizations to cancel. 




loss of concentration around the mean polarizations for the 
two modes is not consistent with the observations, and a 
further assumption is needed to overcome this. Before con- 
sidering how this can be achieved, we relax the assumption 
of orthogonality. 

3.3 Nonorthogonal polarizations 

The example shown in Figure 5 differs from Figure 4 only 
in that the centroid for mode 2 is not orthogonal to that for 
mode 1, and is displaced from the antipodal point by 2/S.il) — 
3°, 2 Ax = 10°. This introduces a favored direction on the 
Poincare sphere, and the polarization points tend to form 
an annulus around this preferred direction. This annulus 
provides a natural explanation for the spreading apparent 
in the observational example in the upper panel of Figure 1. 
However, the absence of a noticeable concentration of points 
around the mean polarization for mode 2 is inconsistent with 
the observations. 

Note that the annulus appears when the ratio of the 
mean intensities is close to unity, but slightly greater than 
unity (fco = 1.1 here), and the variance (A2 = 0.05 here) 
corresponds to a standard deviation roughly equal to the 
difference (A2 ~ (fco — 1)"^). The large spread is due to the 
small variance in this ratio causing there to be a large num- 
ber of weakly polarized bursts due to the mean polarizations 
nearly canceling. 



3.4 Sometime invisible mode 

There is an obvious way of rectifying the absence of con- 
centrations of points around the modal point in Figures 4 
and 5: include a component of the emission that has such 
a concentration. One way in which this can be achieved in 
our model is to introduce one or more additional subsources 
that have such concentrations. However, it is of interest to 
explore an alternative way of achieving the same result with 
only one subsource: assume that some pulses are visible in 
mode 2 but not in mode 1. This is the purpose of our pa- 
rameter pv which specifies the fraction of pulses in which 
mode 1 is invisible. 

In Figure 6 we show how Figure 5 is modified by in- 
troducing the parameter p„ — 0.2, which corresponds to 
mode 1 being invisible in 20% of the pulses. This adds the 
peak around the centroid for mode 2 in Figure 6, compared 
with Figure 5. This peak is due entirely to the pulses for 
which the component in mode 1 is assumed to be invisible. 
In both Figure 6 and Figure 7, there are three peaks in the 
histogram for PA, at -90°, ^ -10° and 45°. The first of 
the three can be seen to wrap around at +90°, due to the 
PA being modulo 180°. The peak at 45° corresponds to the 
cases where mode 1 is invisible. As for the two other peaks 
at negative PA, the picture is similar to what one expects 
for nonorthogonal modes (McKinnon 2003b): they are sep- 
arated by < 90° with a bridge between them. This is due to 
the addition of the non-orthogonal polarizations of modes 1 
and 2, as described for Figure 5. In fact a similar feature 
appears in the histogram for Figure 5, which is not shown. 
Besides the nonorthogonality of the centroids of modes 1 
and 2, the different probability distributions for the polar- 
izations of the two modes is an important ingredient in the 
structure of the PA histogram. 

The annulus first seen on the Poincare sphere of Fig- 
ure 4, remains in Figures 6 and 7. As discussed above, this 
is due to the almost constant ratio of intensities and the dif- 
ferent shapes of the distributions of modes 1 and 2, which 
results in instantaneous nonorthogonal polarizations being 
added. 



3.5 Multiple subsources 

The introduction of the parameter pv has another conse- 
quence that is evident in the fifth histogram in Figure 6: 
there is a narrow peak in the degree of linear polarization, 
corresponding to a small fraction of 100% polarized bursts. 
This is an artefact of allowing a fraction of the pulses to 
be completely polarized in mode 1. This artefact can be 
eliminated by allowing multiple subsources. An example is 
shown in Figure 7, which differs from Figure 6 in that there 
are A'^ = 3 subsources. In this case the parameter p„ = 0.2 
is complemented with pjv = 0.4. 

The distribution of points in Figure 7 is similar to that 
in the upper panel in Figure 1, and we could improve the 
agreement by further adjustment of parameters. However, it 
is also clear that there is no unique way of simulating this 
particular example. We chose to introduce the parameters 
Pv and Pjv to overcome the loss of concentration around the 
centroid of mode 2 when modeling the spread distribution 
for mode 1 with only one subsource. There is a wide de- 
gree of freedom if one assumes that the different subsources 
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Figure 6. As for Figure 5 but with an a fraction pv = 0.2 of the 
pulses having mode 1 invisible. 
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Figure 7. As for Figure 6 but with N = 3 subsources, The 
fraction p„ = 0.2 of pulses with mode 1 invisible does not apply 
to a fraction pjv = 0-4 of the subsources, which have py = 0. 



can have significantly different properties. One could postu- 
late different subsources to produce different features in any 
given observation; specifically, one subsource to produce tfie 
concentration of points for mode 2, with a negligible con- 
tribution from mode 1 (fco S> 1) and a second subsource 
to produce mode 1 spread out around an annulus, with no 
significant focus for mode 2, as in Figure 5. 



4 PHYSICAL INTERPRETATION 

In discussing the physical ingredients in our model we start 
with some general comments, then concentrate on two fea- 
tures of the observations that seem difficult to explain: 
nonorthogonality of the OPMs, and asymmetry between the 
OPMs. We also comment on the difficulty of explaining the 
separation into two modes with comparable intensities and 
the limiting polarization. 



4.1 Physical ingredients 

Important physical ingredients in understanding the polar- 
ization of pulsar radio emission include birefringence, gy- 



rotropy, ray tracing, separation into two natural modes and 
limiting polarization. The first three of these ingredients are 
relatively well understood in principle, although there re- 
main major uncertainties in applying them to pulsars. Bire- 
fringence is present in any magnetized plasma: there are two 
natural wave modes in the medium, and these have different 
refractive indices and different polarizations. Gyrotropy im- 
plies that the natural modes are elliptically (rather than 
linearly) polarized, and in the present context gyrotropy 
is due to differences between the distributions of electrons 
and positrons (Melrose & Luo 2004). Under most condi- 
tions the natural modes are orthogonally polarized, which 
implies that their polarization ellipses are orthogonal with 
opposite handedness. If the medium is inhomogeneous, then 
birefringence implies that the rays corresponding to the two 
natural modes propagate along different ray paths, and ray 
tracing involves determining these paths for specific mod- 
els of the inhomogeneous plasma (Barnard & Arons 1986; 
Petrova 2001; Petrova 2002). 
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4.2 Nonorthogonality and asymmetry between 
OPMs 

Observationally, there are clear examples where the observed 
polarizations are concentrated around two mean polarizar 
tions that are not orthogonal to each other. However, we 

do not regard the nonorthogonality of OPMs as implying 
a nonorthogonality of the natural modes of the medium. 
OPMs refer to the polarizations seen by the observer. Trac- 
ing the rays for each OPM back to the source, the birefrin- 
gence implies that these originated from different points in 
the source, with slightly different angles of emission, such 
that the rays are parallel and coincident at infinity. Physi- 
cally, the nonorthogonality of the OPMs reflects the differ- 
ence in the polarizations of the modes at these points and 
with these directions of emission, and does not reflect an 
actual nonorthogonality of the natural modes at any given 
point in the medium. An estimate of the nonorthogonality 
may then be used to infer a constraint on the properties of 
the emission given a model for the medium through which 
the radiation propagates to the observer. 

Nonorthogonality of the modes is included in two ways 
in the model. Besides the explicit nonorthogonality in the 
centroids of the modes, there is also a nonorthogonality aris- 
ing from the choice of random number from the independent 
distributions of the polarizations for the two modes in each 
simulated pulse. We envisage the source being composed of 
many subsourccs, whose polarizations differ due to small 
differences in relevant parameters between them. Pulse-to- 
pulse variations in both the source regions and along the ray 
paths for the two modes are assumed to lead to a degree of 
randomness in the observed polarization that is described 
by both the centroids for the two modes being different and 
the variations about these centroids being independent for 
the two modes. 

Consider the interpretation of the upper panel in Fig- 
ure 1, which is an example where the two modes are asym- 
metric, with one concentrated around a centroid, and the 
other spread out, with no noticeable concentration around 
the antipodal point. Our model shows that this is consistent 
with two OPMs provided that (i) the OPMs are not strictly 
orthogonal, and (ii) there are two subclasses of pulses with 
(a) the intensities nearly equal in one subclass, leading to 
the spread out distribution, and (b) one mode is invisible in 
the other subclass, leading to the concentrated distribution 
of points. Although some nonorthogonality in the means of 
the modes is required, it is also clear that there are two other 
important ingredients in simulating this case: the spread in 
polarization about the means, and a tight (but not exact) 
correlation of the intensities of the two modes (in the first 
of the two subclasses of pulses). Physically, this emphasizes 
the importance of the pulsar radiation separating into two 
modes that propagate along different rays paths and have 
nearly equal intensities. 

Our interpretation of the annulus in Figure 1 may 
be compared with that suggested by Edwards & Stappers 
(2004), who interpreted it in terms of generalized Faraday 
rotation (GFR). The GFR interpretation requires that po- 
larized radiation be incident on a region in which the polar- 
izations of the natural modes is different from that of the 
incident radiation. Different amounts of GFR in different 
pulses can then lead to a distribution of points around an 



annulus, with the anrmlus oriented orthogonal to the diag- 
onal through the Poincare sphere defined by the polariza- 
tions of the natural modes. In one sense the model based 
on GFR is complementary to the empirical model proposed 
here. The empirical model is based on probability distri- 
butions, and no specific physical assumptions are made in 
writing down the probability distributions. GFR provides a 
physical basis for a probability distribution that, according 
to Edwards & Stappers (2004), reflects the distribution of 
points seen in Figure 1. However, this is not the interpre- 
tation that we suggest. Our interpretation is bsised on the 
polarizations of the two orthogonal modes nearly canceling: 
provided the two modes have nearly eqxial intensities their 
sum is weakly polarized, with a wide spread in polarization. 
With our assumption that the OPMs are 100% polarized, 
if the intensities are markedly different (in a given pulse) 
then the polarization is necessarily close to that of the mode 
with the higher intensity. An implication of the interpreta- 
tion we propose is that the points around the annulus must 
be weakly polarized. There is no such implication with the 
GFR interpretation. 

4.3 Separation into two modes 

There is a major difficulty in understanding how radiation 

can be produced in two modes with nearly equal intensities 
and significantly different ray paths. 

How the radiation becomes separated into two natural 
modes is poorly understood. Most radio emission mechar 
nisms favor radiation into a single natural mode. For exam- 
ple, any maser process that causes one mode to grow faster 
than the other leads, after many growth times, to the faster 
growing mode completely dominating. In principle, this need 
not be the case: if the growth rate is larger than the rate of 
generalized Faraday rotation, and if the maser is intrinsically 
polarized with a polarization different from that of either 
natural mode in the medium, then the growing radiation can 
be an intrinsic mixture of the two natural modes (Melrose 
& Judge 2004). Although the conditions required for this 
to apply seem implausible, the alternatives seem even less 
plausible. The alternative is that the emission mechanism 
results in a single mode, and that the separation into two 
modes occurs somewhere along the propagation path. Al- 
though mode coupling does occur due to inhomogeneity, it 
is usually a weak effect, whereas the interpretation of OPMs 
requires comparable intensities in the two modes. This is es- 
pecially the case for the interpretation of the broad spread in 
polarization shown in Figure 1: our modeling suggests that 
this implies that I2 is strongly correlated with /i. Effective 
mode mixing could occur due to reflection of waves off a 
sharp boundary, which would need to be near the source 
region to be effective. (Far from the source the refractive in- 
dices are very close to unity, precluding significantly different 
ray paths for the two modes, as seems to be essential for the 
interpretation of OPMs.) However, there is no model which 
incorporates reflection off sharp boundaries. Moreover, near 
the source one of the modes should have wave properties 
close to the vacuum, referred to as the X mode by Arons 
& Barnard (1986), and this neither reflects off sharp gra- 
dients nor otherwise couples to the other mode. In brief, it 
is very difficult to account for pulsar radiation that is an 
approximately equal mixture of two natural modes, despite 
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the overwhelming observational evidence that the radiation 
is such a mixture. Hopefully, further use of empirical models 
will lead to information on the polarizations of the OPMs 
that will help constrain the mechanism that leads to the 
separation into two modes. 

4.4 Limiting polarization 

The observed polarization is clearly elliptical in some cases, 
implying that the natural modes of the pulsar plasma are 
elliptical at the point where the radiation effectively escapes 
the magnetosphere. As the radiation in a given mode prop- 
agates through the medium, its polarization adjusts contin- 
uously so that it remains that of the natural mode at every 
point along the ray path. (A small leakage into the other 
mode occurs, implying some mode coupling.) This putative 
point is referred to as the polarization limiting region, which 
may be defined as the region beyond which the medium be- 
comes ineffective in changing the polarization of the radi- 
ation propagating through it. Polarization limiting is most 
likely to occur near the cyclotron resonance where the po- 
larization of the natural modes is changing fastest as a func- 
tion of distance along the ray path (Melrose & Luo 2004). 
Thus the empirical modeling of the polarization provides in- 
formation on the polarization limiting region, and on how 
it varies statistically from pulse to pulse. Such information 
should help to identify the location of the polarization lim- 
iting region. 



5 CONCLUSIONS 

In this paper we present an empirical model that is useful in 
simulating data on the Stokes parameters for observations 
of single pulses at a given pulsar phase. The observed polar- 
ization changes from pulse to pulse, and the polarization in 
any given pulse can be quite different from the (mean) po- 
larization found by summing over a large number of pulses. 
The underlying model for the interpretation of the observed 
polarization is in terms of two OPMs: the observed radi- 
ation is assumed to be a mixture of the two (completely 
polarized) OPMs, with intensities 7i, I2, that are partially 
but not completely correlated. We model this by assuming 
that the intensity for mode 1 has log-normal statistics, with 
the mean intensity set to unity without loss of generality 
and the variance in the natural log set to Ai = 0.8 based on 
observation, (Cairns, Johnston & Das 2004). The ratio, k, of 
the intensities in mode 1 and 2 is assumed to have a gaussian 
distribution with a mean fco and a variance A2. The results 
of the simulations are sensitive to this correlation, and the 
choice of these parameters are severely constrained by the 
data, with fco near but not equal to unity and A2 small but 
nonzero, cf. Table 1. 

In the empirical model presented here, we assume that 
the polarization of each OPM is determined by a prob- 
ability distribution, with a mean and a variance on the 
Poincare sphere. The model allows some non-orthogonality 
in the OPMs. The intensity of mode 1 is assumed to vary 
from pulse to pulse with log-normal statistics. The intensity 
of mode 2 is partially but not totally correlated with the 
intensity of mode 1. When the polarization of each simu- 
lated pulse is plotted on the Poincare sphere, pulses with 



11 ^ I2 and Ji <C I2 lead to peaks around the centroids for 
modes 1 and 2, respectively, and pulses with Ji ~ I2 lead 
to a spreaded out distribution of (weakly polarized) points 
on the Poincare sphere. The model also allows a single pulse 
to consist of emission from a rmmber, N , of subsourccs. In 
the case where we have multiple subsources, we choose the 
statistical parameters (means and variances) of each sub- 
source to be the same. This leads to some difficulties that 
we overcome by a somewhat artificial assumption that allows 
both concentrations around the OPMs and a broad spread 
of weakly polarized points, as required to simulate some ob- 
servations. An alternative (also artificial) assumption would 
be to allow the statistical parameters to be different for dif- 
ferent subsourccs, relying on some subsources to produce 
the broad spread and other subsources to produce the well- 
formed peaks about one or both OPMs. Other correlations 
in the simulations are shown through histograms for the four 
Stokes parameters, and also for the position angle and the 
degrees of linear and circular polarization. Such histograms 
make it possible to describe the variations in the intensity 
and in the degree of polarization, neither of which are rep- 
resented on the Poincare sphere. 

Our model contains a large number of free parameters, 
with this number increasing approximately as the power 
of the number of subsources. However, given the general 
framework of the OPM model, this number of parameters 
is needed to describe a single subsource, and one must rely 
on comparison with observation to reduce or otherwise con- 
strain the number of free parameters. 

The observational examples on whose interpretation 
we concentrate, cf. Figure 1, includes one phase with well- 
defined OPMs and another with a concentration of points 
around one OPM and a diffuse distribution of points that 
extends around an annulus passing only roughly through 
the antipodean point on the Poincare sphere. The spread- 
out case is the more difficult to explain. We find that it 
requires non-orthogonal OPMs with tightly correlated in- 
tensities, with best fits found for the mean and variance of 
the ratio, k = I2/I1, of the intensities having the values 
fco = 1-1, A2 = 0.05. However, this precludes a well-formed 
peak around either mode, and an additional assumption is 
needed to account for both a concentrated peak for one mode 
and a broad spread for the other. We show that this can 
be achieved by assuming that some fraction of the pulses 
contain only one mode. (A physical justification for this as- 
sumption is that the other mode is sometimes refracted into 
a direction that does not intersect the Earth.) The impor- 
tant parameter that distinguishes between the top and bot- 
tom panel in Figure 1 is the variance in the ratio of the 
intensities: for A2 = 0.5 the model gives two well-defined 
OPMs. The explanation is straightforward: a large variance 
in /2//1 implies that most pulses have either I2 <C Ii or 

12 S> /i, giving a polarization distribution close to that for 
mode 1 and mode 2, respectively. This further emphasizes 
the importance of the tight correlation required to account 
for a broad distribution on the Poincare sphere. Note the im- 
plication that the degree of polarization of the points that 
define the broad distribution must be low (very nearly equal 
mixtures of two OPMs). However, we should emphasize that 
although these qualitative conclusions are robust, we do not 
claim that any specific parameter fit is unique. Indeed, it 
seems likely that there is considerable freedom in the inter- 
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play between the variance of the polarization of one OPM 
(012,62) and in the ratio the intensities (fco, A2) in account- 
ing for mode 2 around an annulus, in modeling the lower 
panel in Figure 1 for example. 

Our ultimate objective is to understand the origin of the 
polarization of pulsar radio emission, and hopefully to use it 
to constrain pulsar models. Our results provide strong sup- 
port for a model in which the radiation escapes in two natu- 
ral waves modes that can be significantly elliptically polar- 
ized. Moreover, the intensities of the two modes arc tightly 
correlated in at least some examples that we have consid- 
ered. We note that there is no accepted emission mechanism 
that can produce roughly equal intensities in two orthogonal 
modes, and no known propagation effect that can be effec- 
tive in producing nearly equal intensities in two orthogonal 
modes for radiation initially in one mode. Further ideas on 
how this dilemma might be resolved are needed. 
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